Phase diagram of the Heisenberg spin ladder with ring exchange 
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We investigate the phase diagram of a generalized spin- 1/2 quantum antiferromagnet on a ladder 
with rung, leg, diagonal, and ring-exchange interactions. We consider the exactly soluble models as- 
sociated with the problem, obtain the exact ground states which exist for certain parameter regimes, 
and apply a variety of perturbative techniques in the regime of strong ring-exchange coupling. By 
combining these approaches with considerations related to the discrete Z± symmetry of the model, 
we present the complete phase diagram. 

PACS numbers: 71.10.Hf, 75.10.Jm 



I. INTRODUCTION 

A number of experiments conducted in recent years 
suggest that multiple spin-exchange interactions have 
a significant role in the quantitative description of the 
physics of low-dimensional cuprate compounds, and a 
qualitative one in determining the properties of 3 He 
adsorbed on graphite surfaces. Inelastic neutron scat- 
tering measurements of the spin-wave spectrum in the 
two-dimensional (2D) cuprate system L^CuO^ 1 and in 
the quasi-lD spin-ladder compound La6CasCu2404i, 2 as 
well as two-magnon Raman scattering measurements on 
doped systems related to the latter, 3 indicate the pres- 
ence of contributions from a four-spin cyclic exchange in- 
teraction on the order of 10-20% of the nearest-neighbor 
superexchange. Measurements of the magnetization and 
heat capacity of 3 He films of various fillings adsorbed on 
graphite surfaces 4 ' 5 have been interpreted 6 in terms of 
cyclic three-, four- and higher-spin exchange processes. 

The four-spin exchange term upon which we will fo- 
cus here arises at fourth order in a strong-coupling 
(small t/U) perturbative expansion of the one-band 7 
and three-band 8 Hubbard model in 2D, and has been 
shown in this limit to give the leading correction to 
the nearest-neighbor Heisenberg model. Recent investi- 
gations of the influence of this interaction on physical 
properties have employed perturbative approaches, 9 ' 10 
spin-wave analysis, 11 ' 12 numerical exact diagonalization 
of small clusters, 6 exact diagonalization in combination 
with conformal field theory (CFT), 13 and the density ma- 
trix renormalization group (DMRG) technique. 14-17 The 
majority of these studies is restricted to the regime of 
weak ring-exchange coupling and to systems with only 
a nearest-neighbor Heisenberg superexchange interaction 
in addition to this term. 

Here we consider a generalized model to gain further 
insight into the nature of the phases and phase transi- 
tions within this class of system. We investigate a spin 
ladder which includes antifcrromagnetic Heisenberg leg, 
rung and diagonal, or cross-plaquette next-neighbor, in- 
teractions, as well as a ring-exchange term. This system, 
depicted schematically in Fig. 1, represents the minimal 
model possessing both the possibility of a four-spin cyclic 
exchange interaction and, by virtue of the diagonal cou- 
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FIG. 1: Generalized spin-ladder system with rung, leg, and 
diagonal superexchange interactions, J±, Jm, and J x , and 
four-spin ring-exchange interaction, K. Each site contains a 
spin S — 1/2. 



pling, points in parameter space with non-trivial exact 
solutions. 

The ring-exchange interaction may be considered to 
be composed of different four-spin coupling terms be- 
tween the spins of a plaquette, subject to the special 
condition that the effective leg-leg coupling is equal to 
the rung-rung interaction, and equal in magnitude but 
opposite in sign to the diagonal-diagonal coupling. Re- 
laxing this constraint on the possible four-spin terms re- 
sults in a more general parameter space with still richer 
phase behavior, but we are unaware of studies of this 
model for arbitrary values of all coupling parameters. 
In addition to the ring-exchange term, a further possi- 
ble section of this general parameter space is given by 
the composite-spin representation of a S = 1 bilinear- 
biquadratic chain, 18 which corresponds to the choice of 
equal leg-leg and diagonal-diagonal couplings. For this 
specific choice there is an explicit mapping between the 
spin-1/2 ladder model and the spin-1 chain. In general, 
the different phases of the spin model, and the transi- 
tions between these, may be distinguished by consider- 
ing the expectation values of nonlocal quantities, such 
as the string order parameter first introduced 19 in stud- 
ies of the spin-1 chain and the Lieb-Schulz-Mattis twist 
operator. 20 The differing effects of the four-spin coupling 
terms on string order parameters defined for the spin-1/2 
ladder were analyzed in Ref. 21. 

We begin by considering the points within the pa- 
rameter space for which exactly soluble models exist, 
and then analyze those regimes in the vicinity of these 
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FIG. 2: Phase diagram of the generalized ladder model of Eq. (1), represented as a function of the coupling ratios J±/K, 
Ju/K, and J x /K. The inclined plane indicated by dot-dashed lines is invariant with respect to a Z4 symmetry, and is given 
by the equation J|| — J x = K. The dashed line passing through A and B represents an exactly soluble model with central 
charge c = 1 between A and B, a "rung-singlet" (RS) phase to the left of point A and a ferromagnetic phase (FM) to the right 
of point B. The full line AA' corresponds to a c = 3/2 conformal field theory and belongs to the Babudjian-Takhtajan (BT) 
universality class; it is the exact transition line between the rung-singlet phase and a staggered dimer (SD) phase. The full line 
BB' separates the ferromagnetic phase from a spin-gapped region of the type proposed by Affleck, Kennedy, Lieb, and Tasaki 
(AKLT). The shaded triangle lying below the ^-symmetric plane represents a critical surface which is also of the c = 3/2 BT 
universality class. 



points which are accessible by perturbative techniques 
based on the CFT corresponding to each exact solu- 
tion, i.e. to the regions where the soluble models are 
critical. The ladder Hamiltonian which we study, and 
its associated soluble models, are based on a plaque- 
tte algebra, by which is meant an algebra of operators 
defined on the four-spin plaquette units of the ladder. 
This algebra satisfies quasi-local commutation relations, 
meaning that operators defined on next- (and further- 
neighbor plaquettes commute. Using the commutation 
relations of this algebra we have found exact solutions 
to the Yang-Baxter equation. 22 These solutions provide 
soluble, S£/(2)-invariant, isotropic spin-ladder Hamilto- 
nians, whose excitation spectrum may be either gapless 
or gapped. 22 



Gapless, or critical, solutions arise at points or lines 
of phase transitions, and may be further mapped in the 
continuum limit to CFTs with a central charge c > 
1, or more specifically to Wess-Zumino-Witten (WZW) 
models. 23 Perturbative approaches constructed around 
such WZW models enable one to access regions of param- 
eter space away from these second-order phase-transition 
points. According to the Zamolodchikov c-theorem, 24 the 
addition of perturbative conformally noninvariant terms 
represented by relevant and marginal operators results in 
a flow cither to another CFT with smaller central charge 
c or to a massive phase. A renormalization-group (RG) 
procedure constructed on this basis of these perturba- 
tions may then be used to establish which fixed points 
are stable, and to determine the connectivity, the flow 
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structure and thus the physical properties of different 
regimes of the phase diagram. 

An exact description of the gapped phases is possible 
by using the idea of "words" on the plaquette algebra to 
construct variational ground states. The direct relation 
of this concept to the matrix-product ansatz 25 allows one 
to obtain both the exact ground state in suitable param- 
eter regimes and the lowest excited states, from which 
the lines of phase transitions may also be deduced. For 
parameter regimes in the vicinity of these exact, gapped 
solutions, the absence of phase transitions means that 
many qualitative properties of the associated, gapped 
phases are known. 

Throughout this study we emphasize the role of differ- 
ent symmetries, in particular the symmetry of the 
model, in obtaining a complete understanding of the 
phase diagram. The relevance of this symmetry to the 
ring-exchange problem was first recognized in Ref. 15 in 
the more restricted context of a "duality" transforma- 
tion about a self-dual point in the parameter space of 
the nearest-neighbor ladder model. Here we use the full 
transformation to relate different soluble models, and to 
develop a perturbative bosonization scheme applicable in 
the regime of strong ring-exchange coupling. 

All of these considerations allow us to deduce the com- 
plete phase diagram for the system in the 3D space of 
coupling ratios J±/K, J\\/K, and J x /K, which for the 
orientation of the reader is presented in Fig. 2. The sym- 
bols and lines summarized in the caption of Fig. 2, and 
the phases and transitions they represent, are explained 
in the course of the analyses to follow. Fig. 2 contains 
as a subset [in the line (J±/K = J\\/K, J x = 0)] the 
K > part of the phase diagram obtained in the DMRG 
studies of Lauchli et at, 16 and confirms the presence of 
all of the phases and phase transitions proposed in that 
work. Placing this subset in the context of the full phase 
diagram permits a significant expansion of our under- 
standing of the nature of these phases, and of the first- 
and second-order transitions separating them. Further, it 
provides certain novel, additional phases and transitions 
inaccessible from the restricted parameter space, yielding 
additional insight into the competition of superexchange 
and ring-exchange interactions. 

In Sec. II we introduce the Hamiltonian of the model 
and summarize briefly its associated exact solutions, all 
of which are critical and belong to different universal- 
ity classes. Further exact and variational ground states 
obtained in different regions of parameter space using 
adapted matrix-product wave functions are presented in 
Sec. Ill, accompanied by a discussion of the nature of 
the associated gapped states and their behavior at phase 
boundaries. Sec. IV contains a complete perturbative 
analysis around the four available CFT solutions, with 
emphasis on the properties of the strong-coupling (large- 
K) regime, which establishes the fixed-point structure 
of the phase diagram. Sec. V presents a summary and 
discussion of the results. 



II. LADDER MODELS 
A. Ring Exchange 

The Hamiltonian of the two-leg ladder (Fig. 1) is 
H = J\\ ^(Sj-Sj+i + Tj-Tj+i) 

i 

+ Jx J^(Si-T i+ i + Tj-Sj+i) 

i 

+ J±Y l S i -T i 

i 

+ i^^(P 4 + P 4 - 1 ), (l) 
□ 

where in the last ("ring-exchange") term the four-site 
permutation operator P4 exchanges spins in a cyclic man- 
ner around each elementary plaquette □ of the ladder, 
and is given by the equation 

-P4 + -P 4 1 = Sj-Sj+i + Tj-Tj+i + Sj-Tj 

+ Sj+i-Tj + i + Sj-Tj + i + Tj-Sj+i 
+ 4[(S 4 -S 4+1 )(T,-T 4+1 ) (2) 
+ (Sj-Tj)(S»+i-Tj + i) 
- (Sj-Tj+i)(Sj+i-Tj)] + j. 

We introduce the orthonormal basis 

|0) = ^(|U)-|IT)), |i) = ITT), 

|2> = ^(|U) + IIT», |3) = |||), (3) 

of singlet and triplet states on each rung, and construct a 
set of corresponding projection operators which generate 
the SU(4) algebra 

xf = (\a){0\)u Xfxp = 6^Xr, (4) 
J2 x ? a = 1. [xf,X] 5 ] = {5^Xf-5 a5 xf)5 ij , 

a 

where a, (5 = 0, 1, 2, 3. Spin operators on both chains of 
the ladder may be expressed in terms of this operator 
basis. 22 The ladder Hamiltonian (1) is thus equivalent to 
a generalized, four-state chain with only nearest-neighbor 
interactions. It can be written as 

H = E^ii(^Ui-^+i) + I j x(^ +1 -^ +1 ) 

i 

+ \K{Pl +i - E+ +1 ) (±Jj_ + K)(X™ + Xgj 
+2KXf X^ + \K, (5) 

where P^ i+1 and P* i+1 are respectively permutation op- 
erators corresponding to ladder legs and plaquette diag- 
onals, and the operators E^ i+1 are unnormalized pro- 
jectors on plaquette-singlet states. These operators sat- 
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isfy simple algebraic relations. 26 Their explicit expres- 
sions are 

3 

pi I _ \ " Y a0 Y /3a 
r i,i+l — Z_^i i 
a,/3=0 

= Y,{l-28 a0 ){l-28^)Xfx^, (6) 

a,/3=0 

and 

^M+i=2|^ +1 )(^ +1 |, (7) 

where the functions are plaquette-singlet states, 

i.e. states with total spin per plaquctte S* 2 = (S x + S 2 + 
S3 + S4) 2 = 0, and may be expressed in terms of the rung 
states \a)i and as 

l^i+i) = \{\0) l \0) l+1 ±\2) l \2) l+1 

T|l>i|3) i+1 T|3)<|l) i+ i}. (8) 

There are precisely two possible plaquette-singlet states, 
corresponding respectively to Ef i+1 and B t ~ i+1 . 

B. Integrable Models 

The complete set of exactly soluble isotropic ladder 
models with short-range interactions has been found in 
Ref. 22 using the algebraic Bethe-Ansatz method (for 
a review see Ref. 30). Not all of these models are in 
the same class as that specified by Eq. (1), and here 
we present those cases relevant to a system with ring- 
exchange interactions. All of these are critical, and 
may be distinguished according to their central charge 
c, which characterizes the different universality classes 
close to integrable points. 

1. c = 2 

The exactly soluble models of this type have the simple 
form 

H = Y J {P l -E i ). (9) 

i 

where Pi and Ei denote respectively p}\ +1 , P* i+ i and 
We thus obtain four different soluble 
Hamiltonians which correspond to the four combinations 
{PV,E+}, {Pl l ,Er}, {P*,E+}, and {P?,Er}. Two of 
these are trivially soluble, because J2i(Pi ~ ^i) 1S tnc 
Hamiltonian of two decoupled chains and $^(-Fj X — E^) 
is the same pair of chains intertwined by the transforma- 
tion S 2i <-> T 2 i,S 2 i+i <-> S 2l+1 ,T 2i+1 <-> T 2l+ i [Eq. (5)]. 
The Hamiltonian of two decoupled chains is represented 
on the phase diagram of Fig. 2 by the point J\\/K = 00, 



with all other couplings equal to zero. Nonzero values of 
these couplings generally induce a spin gap, but because 
of the competing nature of different interaction terms 
there is a possibility of critical behavior for some combi- 
nations of these couplings. This will be demonstrated ex- 
plicitly in Sec. IV. Of the remaining two solutions within 
the general ansatz (9), one is given by the combination 
{P i , which in terms of spin operators is 

Hd = SjSj + i + T/Tj_|_i + S/Tj + i + TjSi + i 

i 

+4[(S i S i+ i)(T i T i+ i) - (S i T i+ i)(S <+r T i )]. 

(10) 

The last, {P* , Ef} is obtained by the same intertwining 
transformation from the previous Hamiltonian, and takes 
the form 

Hdi = ^ Sj-Tj+i + TiS i+ i + S^Sj+i + TYTj + i 

i 

+4[(S i T i+ i)(T i S i+ i) - (S^ +1 )(T,T i+1 )]. 

(11) 

For all four masslcss models one obtains two decou- 
pled Bethe-Ansatz equations which correspond to the 
D 2 = SU(2) x SU(2) algebra. The energy of the Hamil- 
tonian is the sum of two expressions for each SU{2) 
component, 31 and thus the Hamiltonian (10) is equiv- 
alent to two decoupled chains. This model is critical 
(no spin gap), with a conformal charge c = 2. Eq. (10) 
constitutes a part of the plaquette term, Eq. (2), while 
we will demonstrate explicitly below that the remain- 
ing contributions to P 4 + P^ 1 induce a spin gap, and 
therefore that the origin of coordinates in the phase di- 
agram of Fig. 2 represents a gapped system. However, 
as shown in Sec. IV, the competition between different 
interaction terms of the Hamiltonian (1) drives the sys- 
tem to a second-order phase transition described by a 
CFT with c = 3/2. The approximate form of this critical 
surface, represented in Fig. 2 by the shaded triangle, is 
obtained in Sec. IV. 

Exact solutions in this class then provide two massless 
models, one in the limit of weak K, and one in the strong- 
K limit, both of which correspond in the continuum limit 
to c = 2 WZW models. The four solutions of the (9) 
are related by two transformations, one of which is the 
intertwining transformation shown above. The other is 
less transparent in nature, and is found by first noting 
that the two unnormalized projectors E~l and E~ are 
related by the X-operator transformation 

X° k a -+ -iXl\ Xf - l Xf, (12) 

for a — 1,2,3 and k = 1,...,N, where N is the total 
number of sites. This is a unitary transformation gener- 
ated by the operator 

U(n/2)=eM-^f2( X °k% (13) 
fe=i 
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which has the property that 

H D = UH 2c U\ (14) 

where H 2c is the Hamiltonian of two uncoupled spin-1/2 
Heisenberg chains. The nature of this transformation 
may be understood from the observation that it maps 
one of the plaquette-singlet state defined in Eq. (8) into 
the other, and thus that it transforms the projectors (7) 
according to <-> . This transformation has the 
additional property that [U(ir/2)] 4 = 1, and therefore 
U is one of the generators of the Z 4 transformation as- 
sociated with the center of 5C/(4). In the next section 
we will demonstrate that this symmetry is essential for a 
complete understanding of the phase diagram. 

2. c = 1 

A second exact solution 22 has the same eigenspectrum, 
but not the same degeneracies, as the spin-1/2 Heisen- 
berg chain. The corresponding spin Hamiltonian, 

H 1/2 = -Kj2(S i -T i+1 +T i -S i+1 ) (15) 

i 

+ ±^£(P 4 + Pf 1 ) + Jx^^Si-Ti), 

□ i 

commutes with the generators ^ - Xf b for any a, b = 
1,2,3, which form an SU (3) subalgebra within the SU (4) 
algebra generated by the A-operators. In view of com- 
pleteness relation in Eq. (4), H X j 2 commutes with the 
total number operator J2i f° r singlets, and any mul- 
tiple of this term may be added to the Hamiltonian with- 
out spoiling the integrability. The full symmetry of the 
model is thus SU(3) x U(l) = U(3). In the framework of 
the equivalent Heisenberg chain, inspired by the Bethc- 
Ansatz solution, the model may be considered as a chain 
of effective SU{2) "spins" L, with 

L\ = i - X, 00 = i - S r TV (16) 

The term J±/2K J2i Sj-Tj in H 1 / 2 therefore corresponds 
to a coupling with a magnetic field h = J±/2K. 32 It is 
well known that the spin-1/2 chain in a magnetic field 
develops an incommensurate critical phase for \h\ < 2 
and has a massive phase for \h\ > 2 (see for example 
Rcf. 33). 

One observes from Eq. (15) that this integrable model 
corresponds to the line in parameter space Jm = 0, 
J x = -A, for arbitrary Jj_. When (\J±/K\) < 4 the 
model is in the critical, incommensurate phase, whereas 
for J±/K > 4 the model has a gapped, rung-singlet 
phase, and in the region J±/K < —4 it has a gapped, 
rung-triplet ground state. The nature of the gapped 
phases and the physical properties of the incommensu- 
rate phase will be discussed in more detail in Sec. III. 
In the critical region the model is described by a CFT 
with central charge c = 1. However, the presence of 



conserved charges generating the U(3) symmetry has the 
consequence 3 that the model possesses additional zero 
modes, and therefore while the conformal dimensions are 
those of the c = 1 theory for the Coulomb gas, the de- 
generacies are altered accordingly. The critical region is 
represented on the phase diagram of Fig. 2 by the thick, 
dashed line AB, while the rung-singlet and ferromagnetic 
rung-triplet regions are denoted respectively by RS and 
FM. In Sec. IV (Fig. 10) we will demonstrate that there 
exists a finite critical region in the vicinity of the critical 
line. 



3. c = 3 

A further exact solution exists which corresponds sim- 
ply to local Hamiltonians proportional to permutation 
operators Pj' and P* . The Hamiltonian corresponding 
toff, 

H\\ = ^2 + Tj-Tj + i + 4(Sj-Si + i)(Tj-Tj + i) 

(17) 

is invariant with respect to the full SU(4:) group and 
describes a soluble "spin-orbital" model. 35 In the con- 
tinuum limit, Eq. (17) corresponds to an SU(4) WZW 
model at level k — l, 36 ' 37 and therefore has central charge 
c = 3. This Hamiltonian commutes with the operator 
^(Si-Ti), and thus is integrable for arbitrary values of 
Jx- One may also add a rung-rung coupling interaction 
to obtain the Hamiltonian 

H{ { =H n + Jx^S i -T i + 2*T^(S r T i )(S i+r T i+1 ). 

i i 

(18) 

which was considered in Ref. 38 for arbitrary J± and for 
two values of K (K = and K = — 1). It was shown that 
for both values of K, three phases appear as a function 
of Ji, and are separated by two quantum critical points, 
3 C + and J c _. For K = and Jl > J% = 2 the model is 
in the rung-dimerized phase, for Jl < J < there is a 
gapless phase with c = 3, and for J < Jt_ ~ —1.79 there 
is a gapless phase with c = 2. For K = — 1 the critical 
points occur at 38 = 1/2 and Jt_ ss —1.29. Quite gen- 
erally, the extent of the c = 3 critical region is expected 
to be strongly reduced by a negative rung-rung coupling, 
while it will expand for positive K. The critical behavior 
of is thus expected to persist for finite regions of pa- 
rameter space, and the position of the critical points to 
depend strongly on the rung-rung coupling. We empha- 
size that the model of Eq. (18) does not appear in the 
phase diagram of the ladder model under consideration 
(1) but we have introduced its properties here for use in 
Sec. IV as one appropriate basis model for a perturbative 
expansion. Similar considerations apply to H x , which is 
obtained from Hn by using the intertwining transforma- 
tion of Sec. IIB.l. 
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FIG. 3: Order parameters of the different phases related by 
the Za transformation within the plane J\\ — J x = K. 



III. EXACT GROUND STATES 
A. Z A Plane 

In this section wc apply certain discrete symmetry con- 
siderations to obtain further insight into the structure 
of the phase diagram. We employ the matrix-product 
ansatz to find exact ground states, and also the bound- 
aries between these, even in the region of parameter space 
where the Hamiltonian (1) is not exactly integrablc. 

We begin by noting that for arbitrary Jj_, when 



J\\ - J-x 



K 



(19) 



the Hamiltonian (1) is invariant with respect to the Z^ 
transformation (13) which generates the mapping Ef — > 
Ef . The condition (19) defines a plane in the 3D space 
of the phase diagram, marked in Fig. 2 by parallel dot- 
dashed lines. The Hamiltonian of the system in the Z^ 
symmetric plane commutes with the singlet total number 
operator J2i 

This Z^ transformation is precisely the one noted 
in Refs. 15 and 39, in which it was referred to as a 
duality transformation. It is a canonical transforma- 
tion for the lowest (s=l/2 x ,s=l/2 ) representation of 
SU(2) x SU(2), and conserves the values of the Casimir 
operators in this representation. In terms of the original 
spin variables it takes the form 



Si — 

Ti = n(Sj 



Tj) — Si x Tj, 
Tj + S^T,. 



(20) 



When expressed using the variables Si,Ti, the Hamil- 
tonian Hd (10) is exactly that of two decoupled spin 
chains. It is clear from the definition (2) of the plaquette 
term that Hp constitutes a part of the operator P^+P^ 1 , 
while the remaining contributions to this term couple the 
two effective spin-1/2 chains to produce a system with a 
finite spin gap. 

In the Z± plane this transformation acts as a sym- 
metry between different order parameters. It maps the 
order parameter (Si — Ti) for antiferromagnetic rung 




FIG. 4: Representation of the rung-singlet (RS) phase. The 
ellipses denote a singlet state of the spins on each rung. 



spin correlations into the vector-chirality order parameter 
(SjXTj) = (Si — Ti), while the leg dimer order parameter 
(Sj-Sj+i — Tj-Tj+i) is mapped into the order parameter 
((Si + TiX-Si+i x T j+ i) + (Si x Ti)(-S i+ i + T i+1 )) for 
scalar chirality (Fig. 3). Phases characterized by order 
parameters related under the Z^ transformation within 
each pair are therefore symmetric within the plane de- 
fined by Eq. (19). The point J± = Jn = K considered in 
the DMRG studies of Refs. 15 and 16 lies in this plane 
(the star in Fig. 2), giving the properties of duality ob- 
served in the results of both analyses. 

From this consideration we may deduce the nature 
of a candidate ground state as a product of linear su- 
perpositions of the two plaquette-singlet states. The 
complete set of plaquette states with fixed angular mo- 
mentum contains one quintuplet (j = 2), three triplets 
(j = 1) and two singlets (j = 0). The two plaquette- 
singlet states are those constructed above, the states 
with j = 1 are created by applying a linear combina- 
tion of operators P-] i+1 , P£ i+ i, and X^Xf^, and the 
state with j = 2 is created by the action of a combina- 
tion of the operators Pf i+ \, Pfi + \, an d Ef i+1 . The ex- 
pressions for these states allow their identification with 
matrix-product-ansatz states, while the explicit form of 
the Hamiltonian (5) corresponds to the operator basis 
for matrix-product-ansatz Hamiltonians. 40 One observes 
that in the ^-symmetric plane, J\\ — J x = K , the con- 
dition (J±/K) > 4 defines a region in the phase diagram 
(Fig. 2) with an exact rung-singlet (RS) ground state, 



Mrs) = |0)i|0) 2 ...|0) JV _i|0) JV , 



(21) 



which is represented schematically in Fig. 4. On the line 
J x = — K , J|| = 0, which lies in the ^-symmetric plane, 
this statement is supported by the exact solution of the 
Hamiltonian i?i/2 (15) of the previous section. From the 
results of Ref. 40, the line J±/K = 4 in the ^-symmetric 
plane is a line of second-order phase transitions into a 
spontaneously dimerized phase of staggered leg dimcr- 
ization (SD). 




FIG. 5: Typical ferromagnetic (FM) ground state for large, 
negative J±. 
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FIG. 6: Staggered leg-dimerized state (SD). 



Similarly, the condition J±/K < —4 defines an exact 
ferromagnetic rung ground state (FM), 



IVfm) = XfXf...X%>_ 1 X%>\0), 



(22) 



in which |0) denotes the global state of singlets on every 
rung. This state is represented in Fig. 5. The tran- 
sition on the line J±/K = —4 is of first order. Both 
transition lines are determined exactly from the points 
where the dispersion relations of the elementary excita- 
tions become massless. In the rung-singlet phase this 
excitation is a propagating rung triplet, while in the fer- 
romagnetic phase it is a rung singlet. The second-order 
transition from the rung singlet phase to the staggered 
dimer phase (Fig. 6) lies in the universality class of the 
spin-1 bilinear- biquadratic chain, which in the continuum 
limit is described by a c = 3/2 CFT, 23 and is a transi- 
tion of the Babudjian-Takhtajan type. This transition is 
characterized by a spontaneous breaking of discrete Zi 
symmetry, which in the ladder is related to translation 
by one a lattice unit. 

From the Z4 duality of the system we deduce that 
at large K (near the origin of Fig. 2) there is another 
second-order phase transition of the same universality 
class from a scalar-chirality phase (SC) (Fig. 7), which is 
dual to the staggered dimer phase, to a vector-chirality 
phase (VC) (Fig. 8) which is dual to the rung-singlet 
phase. In the next section we provide further arguments 
in support of this statement. Figs. 7 and 8 show respec- 
tively typical static configurations of spins in scalar- and 
vector-chirality phases; these should be understood only 
as indicating the preferred instantaneous spin configura- 
tions, which in fact fluctuate rapidly, such that the phases 
exist only in the sense of finite average values of the cor- 
responding spin correlation functions. All four phases 
have only short-ranged correlations in the chain direc- 
tion. These states illustrate the increasing dominance 
of the K term, which favors configurations in which all 
spins on a plaquette are mutually perpendicular to maxi- 
mize their solid angle, 11 in the competition with nearest- 





FIG. 8: Schematic representation of a static vector-chirality 
configuration (VC). 



neighbor exchange terms which favor antiparallcl spins. 

Finally, the first-order phase transition from the fer- 
romagnetic phase takes the system to a type of Afflcck- 
Kennedy-Lieb-Tasaki (AKLT) state 41 formed by effective 
spin-1 variables represented by the rung triplets (Fig. 9). 
It differs from the ferromagnetic phase in that coherence 
between the individual rung triplets is not established 
until the transition at J±/K = —4. This state is ob- 
tained directly from the vector-chirality phase in a first- 
order transition at small, negative J± . We stress that the 
letters indicating the locations of these phases in Fig. 2 
should be understood to refer only to the plane of exact 
Z 4 symmetry including the rung-singlet and ferromag- 
netic phases which lie outside the region marked by the 
dot-dashed lines. However, because there are no phase 
transitions other than those present on the diagram, the 
gapped phases obtained for model parameters outside 
the plane are connected continuously to those of the ex- 
actly known states, and their physical properties evolve 
continuously as the values of the interactions are further 
changed away from the Z4 plane. 



B. Incommensurate Line 

With these results, a heuristic understanding of the 
origin of the incommensurate phase appearing on the 
critical line AB in Fig. 2 (Sec. II) may be obtained by 
considering the role of J± as an effective magnetic field. 
We stress first that the incommensurate phase arises with 
no breaking of 51/(2) symmetry in the space of the spin 
variables Sj and Tj. The variable conjugate to the ef- 
fective field, ^ Sj-Tj, corresponds [Eq. (16)] to an aver- 
age singlet density which varies from at the boundary 
to the ferromagnetic phase (Jj_ — > —4K) to 1 at the 
rung-singlet phase boundary (Jj_ — > AK), in exact anal- 
ogy to J2i $i f° r the spin-chain problem. For the lad- 
der, Sj-Tj also expresses the spin correlation on each 
rung, which may be characterized by defining a variable 



FIG. 7: Representation of a spin configuration with finite 
scalar-chirality (SC) correlation function. 



FIG. 9: AKLT-type ground state. The dashed ellipses repre- 
sent effective spin-1 variables formed by triplet states on the 
ladder rungs. 
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k± = cos _1 ((S i -T i )/|S i ||T i |). This effective wave vector 
across the ladder varies continuously from to tt over the 
range —AK < J± < AK, and takes the value k± = tt/2 at 
J± = 0, where the ring-exchange term, favoring locally 
perpendicular spins, is dominant. 

The incommensurate state in the ladder direction is 
crucially dependent on the parameters of the c = 1 
line, notably the special role of the J x term. Only for 
J x = — K is the spectrum massless, with quasi-long- 
ranged correlation functions classifiable by a ladder wave 
vector q which reflect the competition between nearest- 
neighbor exchange and ring exchange, again in exact 
analogy with the spin-1/2 Hcisenberg chain in a mag- 
netic field. 42 The nature of the incommensurate correla- 
tions arising from the excitation spectrum is elucidated 
by considering the effective spin variables Lj. The anal- 
ogous "spin" correlation function is given by 

(Li-Lj) - (Xf a Xf + XfX° a ) 

+2<(I-X°°)(I-X° )>, (23) 

where the first line corresponds to terms of the form 
LfLJ + LJ L~t and the second to L*L Z 4 . In terms of 
the original spin variables one obtains 

(Ivly) = |<(S-T) i (S-T) J -) + 2<(SxT) i (SxT) j ) 
+ i((i+2S r T i )(i+2S j .T j )). (24) 

Thus the dynamical correlation function exhibiting a 
peak at the incommensurate wave vector q is in fact the 
sum of the staggered rung, vector chirality, and singlet 
density correlation functions, although from the massless 
nature of the excitations at q one may expect peaks in 
each function individually. Again the terms in the first 
line correspond to L^LJ + L^Lj~, and are symmetrical 
under the Z4 transformation (Fig. 3) in the plane: their 
peaks appear at the same wave vector for all values of Jj_ , 
but their intensities may differ. The effective LfLj terms 
(second line) are not symmetrical with the others for all 
coupings Jj_ 7^ (finite effective field) , but as in the spin 
chain may be expected to show the same continuous evo- 
lution of the incommensurate peak position from q = 
at J± = — AK to q = it at J± = 0, and back to q = at 
J± = AK. We stress again that this behavior is a spe- 
cific property of the incommensurate line and cannot in 
general be expected to be clearly visible in other parts of 
the phase diagram not in its immediate vicinity, including 
the regions corresponding to DMRG analyses performed 
to date. However, in Sec. IV we isolate a region in which 
the incommensurate properties (albeit for a gapped exci- 
tation spectrum) may persist over a significant range of 
parameters. 

IV. RENORMALIZATION-GROUP ANALYSIS 

In this section we study the perturbations around dif- 
ferent CFT solutions, which either appear in the phase 



diagram of the model (1) or are closely related to it. We 
adopt a variety of techniques to address the nature of the 
ground states in regions away from these soluble points, 
and to establish the fixed-point structure of the phase di- 
agram. Because all of the transition points and lines are 
known, the relevance of the operators within Eq. (1) as 
perturbations of the exact solutions determines the flow 
under renormalization and thus the dominant physical 
properties of the intermediate regions. In each subsec- 
tion and in the Appendix, known results are summarized 
briefly while those which are new in the current context 
are presented in detail. 

We have found four models suitable for this type of 
analysis. From the exact solutions of Sec. II, there are 
two CFTs with c = 2 which correspond to two decoupled 
chains described by the soluble Hamiltonian of Eq. (10); 
one solution is relevant for weak K, and the other for 
strong K. Here we concentrate primarily on the strong- 
coupling regime, i.e. large K, for which significantly less 
is known: the fermionized version of the limiting model of 
two effectively decoupled spin chains is analyzed in sub- 
section A by a conventional RG procedure. We have ar- 
gued from the presence of the Z4 symmetry (Sees. II, III) 
that in the strong-coupling limit there is a second-order 
phase transition with c = 3/2, represented by the shaded 
triangle on the phase diagram of Fig. 2. The Zamolod- 
chikov c-theorem 24 then demands that weak- and strong- 
coupling limits are disjoint in RG sense, i.e. that there is 
no continuous flow from weak to strong coupling. To de- 
scribe the intermediate regime between the two second- 
order phase transitions with c = 3/2 we therefore use 
in subsection B another CFT with c = 3, which is de- 
scribed by the model (18). By perturbative analysis of a 
fermionized model we find that this intermediate region 
is gapped; an alternative perturbative treatment of the 
SU(A) model (17) has recently provided similar results. 39 
Finally, the same approach may also be applied to the 
c = 1 CFT which corresponds to the solution (15). In 
the vicinity of the line AB on the phase diagram of Fig. 2, 
we employ rather general arguments in subsection C to 
reveal the presence of a massless region in one plane of 
the phase diagram. The CFTs with c = 3/2 thus rep- 
resent unstable points in parameter space and must be 
accompanied by a flow towards stable fixed points. Based 
on the arguments of Ref. 23 the stable theory should be 
a SU{2)k=i CFT with c = 1, the natural candidate for 
which is the CFT corresponding to the model of Eq. (15). 



A. c = 2 CFTs 

A weak-coupling bosonization analysis is appropriate 
for the limit of two quasi-decoupled chains when Ji | 3> K 
and J11 3> Jj_. The ladder system with biquadratic ex- 
change has been shown to undergo a second-order phase 
transition at which the behavior of the massless modes is 
governed by a CFT with central charge c = 3/2. 43 The 
consideration of a ring-exchange term is technically iden- 
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tical (see Ref. 10), and the results are found to be in good 
agreement with those from other approaches. 

When performing a perturbative analysis around 
CFTs, all contributions to the spin Hamiltonian should 
be classified according to the scaling dimensions of the 
operator content of the corresponding CFT. In the case of 
perturbations around the limit of two decoupled chains, 
in both weak- and strong- if limits, four-spin interaction 
terms and interleg couplings are expressed in terms of two 
fundamental c = 1 WZW fields, g a and g' a (a=0,l,2,3), 
with conformal dimension and of Kac-Moody cur- 

rents J and J with dimensions (1,0) and (0,1). This 
theory can be further expressed in terms of four differ- 
ent Ising models, 43,44 i.e. of order-disorder fields, energy 
operators and Majorana fermions. The representation of 
these fields by four Ising models is summarized in the 
Appendix. The connection to the spin variables is given 
by 



Si 
S(x) 
T(x) 



aS(x) 

3 + 3 + {-l) x/a eTi[g(x)a} (25) 
J'+J' + (-lf /a eTr[</(x)cr]. 



We note that in Eq. (25) the uniform and staggered parts 
of the spin-density operators have different conformal di- 
mensions, and that 9 is a nonuniversal normalization 
constant. 

The expression in terms of Ising-modcl fermions for 
the general ladder Hamiltonian, which includes arbitrary 
leg-leg, diagonal-diagonal and rung-rung couplings, is de- 
rived in the Appendix, and is valid around the limits 
of two decoupled chains. There are two such limits, 
the weak- if regime which corresponds to the two initial 
chains, and the strong-if regime which corresponds to 
two ^-rotated chains (14), resulting in the Hamiltonian 
(10). The expressions of the Appendix are valid for both 
cases. 

The resulting continuum-limit Hamiltonian for the sys- 
tem in weak- and strong-coupling regimes may be ex- 
pressed in terms of four Majorana fermions with different 
singlet and triplet masses, 43 



H = 



+ v t {rd x r-rd x r)+m t rr\ 

a=l,2,3 

+ffmarg. (26) 



The Hamiltonian of the marginal interactions, 



,„ M ,. g - . dx[\iOi + A 2 2 ], 

Ox = vW 2 ^ 2 + V^V 3 ^ 3 + V' 2 V ;2 V' 3 V' 3 , 

02 = V-VW^ 1 +V> 2 ^ 2 +V> 3 ^ 3 ), (27) 



contains current-current contributions of the forms (Ji + 
Ji)( - J2 + J2), arising from the interleg interactions, and 



— (J1J1 + J2J2) from the intraleg couplings, as well as 
a contribution from the normal-ordered marginal prod- 
uct Tr(a a g)Tr(a a g')Tr(a b g)Tr((T h g') which originates 
in the four-spin term (Sj-Tj)(Sj-|_i-Tj_|_i). Explicit ex- 
pressions for the marginal couplings Ai and A 2 are given 
the Appendix. 



Weak Coupling 

In the weak-coupling regime the results of the Ap- 
pendix provide the expressions 

v2 , 



m t = J±- 2 J x - 20\ £ K 
m s = — 3 J_i_ + 6 J : 



12A 2 if, 



(28) 



for triplet and singlet masses, where A is another nonuni- 
versal quantity. The marginal interactions rcnormalizc 
these masses, and the phase transition occurs when the 
renormalized triplet mass vanishes, i.e. m\ cn = 0. On the 
other hand, the matrix-product ansatz (Sec. Ill) gives the 
phase transition line exactly. One may therefore attempt 
to specify the rcnormalization by taking (for example at 
J x = 0) the value of A 2 to be in agreement with the exact 
result, which suggests that A 2 = 1/5 for this second-order 
phase transition. 



Strong Coupling 

In the strong-coupling regime we analyze perturbations 
around the CFT which corresponds to the exact solution 
of Eq. (10). For this it is convenient to perform the trans- 
formation (20), which from Eq. (1) yields 

H = ^K(^-S J +i+T r f m ) + (J ± + if)(S i -f i ) 

i 

+\Jx {Si*T^ + i+Ti-Si+i + Sj-S^+i+Ti-Tj_|_i 
-4[(S i -S i+1 )(T r f - (S r T <+1 )(S i+r T i)}} 
+ 2^||{S»"Sj+i+Ti-Ti+i+Sj-Ti+i+TVSi+i 

+4[(Si-S i+ i)(T i-f <+ i) - (Si-f i+ i)(S i+ i-T ;)]} 
+2i^(S r T0(S» + rT»+i). 

i 

By substituting 



(29) 



in (A9) the expressions for triplet and singlet masses be- 
come 

mt = J_L+K-J x -J n - A 2 (16J|, + AK - 16J X ), 
m s = -3 J ± - 3if + 3 J x +3Jm +12A 2 if. (30) 



Jl 


= K/2, 


Jd 


= J x /2 + J,|/2 


Vll 


= ~2J X + 2J||, 


Vdd 


= 2J X - 2J||, 


Vrr 


= 2K, 


Jr 


= J±+K, 



10 



The marginal current-current interactions act again to 
renormalize these masses. Setting J x = and taking 
J± = J|| = J yields the value of the second-order phase 
transition point as K/J — 16A 2 /(1 — 4A 2 ). A comparison 
with the DMRG results of Ref. 16, which suggest that at 
strong coupling the transition point is Kj J w 5 — 6 [in 
the units of Eq. (1)], yields A 2 w 1/7, a value rather close 
to that obtained in the weak-coupling regime. After this 
the equation mt = gives the form of the critical surface 
in the strong-coupling limit, which is represented by the 
shaded triangle in the phase diagram of Fig. 2 and by 
the bold, dashed line in Fig. 10; the shaded triangle lies 
below the ^-symmetric plane, but does not include the 
origin of the coordinate systems. We note here that the 
perturbative considerations applied above are in a strict 
sense questionable, because in the problem under consid- 
eration the "perturbative" rung-rung coupling is of the 
same order as the coupling corresponding to the unper- 
turbed CFT, which is K/2. For a ladder system in which 
Vrr is not fixed to be 2K, and is small in comparison 
with Vll and Vdd, the second-order phase transition 
with central charge c = 3/2 in the strong-coupling limit 
is accessible by perturbative analysis around the exact 
solution (10). Thus for consistency this term should be 
added to the Hamiltonian with a small coupling constant, 
Vrr, in order to study the RG behavior. We remark also 
that for J x = Jii = J, the leg- leg and diagonal-diagonal 
biquadratic terms in Eq. (29) cancel explicitly and the 
effective theory is described by massive fermions with a 
singlet-triplet mass splitting generated by Vrr and Jr. 

The lines of phase transitions correspond to the vanish- 
ing of the triplet mass m t ■ The theory is then equivalent 
to three massless fermions and is therefore described by 
a c = 3/2 WZW model. Because of the SU(2) symmetry 
of the ladder model (1), this is a SU(2) WZW model at 
level k = 2. The bare singlet and triplet masses acquire a 
rcnormalization due to marginal interactions which can 
in principle be computed from the RG analysis. By con- 
sidering the operator product expansion (OPE) for oper- 
ators Oi, O2, VWo, and -0a (a = 1, 2, 3) (see Appendix 
and Eqs. (25,27)), we deduce the one-loop RG equations 



d{\nL) 

dg s 
dQnL) 

dg t 
dQnL) 



2tt(A 2 + A 2 ), 



d\ 2 



d(lnL) 



47tAiA 2 , 



g s + 3irg t \ 2 , 

g t (l + 2n\ 1 ) + irg s \ 2 



(31) 



for the marginal couplings Ai and A2 , and for the singlet 
and triplet couplings g s and gt, which are respectively 
proportional to — m s and mt- A similar analysis is per- 
formed in Rcfs. 43 and 44). The first and the second 
equations are decoupled from the others, and can be in- 
tegrated (in variables A± = ± A 2 )) to yield 44 



A± = 



A 



(0) 



From this it follows that for negative initial A± these cou- 
plings are marginally irrelevant (and renormalize to 0), 
whereas for positive values they are marginally relevant 
(and develop exponential gaps). The initial (ultravio- 
let) fixed point is given by g s = g t — 0, which corre- 
sponds to the CFT with c = 2 (two uncoupled chains), 
and the RG equations for g s and g t yield the renormal- 
ization of masses under the RG flow. The new unstable 
fixed point is then defined from the renormalized value 
of m t as m* t = 0. Because of these renormalization ef- 
fects the exact transition line is difficult to estimate, but 
the presence of this transition is to be expected from Z4 
symmetry arguments. 



Self-dual model 

From the analysis above one may conclude that the 
lines (or surfaces) of the two second-order phase transi- 
tions with c = 3/2 are not ^-symmetric, because the 
rung- rung coupling is proportional to K. The Hamil- 
tonian which would have perfect symmetry of the two 
second-order phase transitions is 

H = Jii^Si-Si+i+Tj-Tj+O + Jfl^Si-fj), 

i i 

+|A-^(S i -S i+ i+T i -f i+ i) (33) 



+V Rfl ^(S r f i )(S i+r f i+1 ) 



1 -87rA^ 0) In(L/L )' 



(32) 



which corresponds to the model with Jr and Vrr inde- 
pendent of K. In the small- Ji | limit we have obtained two 
SU(2) WZW models perturbed by relevant and marginal 
interactions which have their origin in the K term and in 
the small couplings Jr, and Vrr. In the strong- K regime 
there are again two SU{2) WZW models perturbed by 
the small Jm, Jr, and Vrr terms. We note that this 
Hamiltonian may be considered as a two-chain Hamil- 
tonian Y^i Pi ~ &i perturbed by terms proportional to 
p\ ' — , or as a ^-transformed Hamiltonian with the 
^-transformed perturbation. Qualitatively, the result- 
ing behavior is the same in both cases. In the continuum 
limit this perturbation is represented by two four-fermion 
terms [see Eq. (41)], one of which gives the marginal con- 
tribution while the other gives the relevant contribution. 



Intermediate coupling 

From the Zamolodchikov c-theorem, 24 a CFT per- 
turbed by relevant and marginally relevant interactions 
will flow cither to another CFT with smaller central 
charge or to a massive phase. For small K and J±, per- 
turbation around the limit of decoupled chains generates 
a flow to a c = 3/2 CFT which is in the universality class 
of the bilinear-biquadratic 5=1 spin chain, a result 
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established recently in Rcfs. 10 and 13. In the strong- 
coupling limit (large K) one may consider the J|| term 
as a small perturbation and the same arguments are ap- 
plicable (but now for the variables S and T), the corre- 
sponding RG flow being from a critical c = 2 CFT to a 
c = 3/2 CFT. The universality class of this transition is 
the same as above: it is described by an SU (2)-symmetric 
CFT with level- fc = 2 Kac-Moody algebra. 

From this one may conclude that the weak-coupling 
regime is not related continuously to the strong-coupling 
regime, and that there is a crossover between the two. 
The natural candidate for this crossover region would 
be the ^-symmetric plane, J|| — J x = K, probably 
in the vicinity of the point J\\ = Jj_ = K on the line 
J x = 0. However, it is important to address the question 
of whether this crossover is a phase transition or a con- 
tinuous change. As shown in Ref. 23, the k = 2 SU{2) 
WZW model is unstable in the sense that it contains rel- 
evant operators which induce a flow to a stable k = 1 
SU(2) WZW, and one would therefore expect another 
second-order phase transition in the universality class of 
the c = 1 WZW model. The alternatives to this sce- 
nario are a first-order phase transition or a continuous 
crossover; a definitive statement is not possible on the 
basis of the present considerations alone. 

Because this special point or line cannot be expected to 
be accessible by RG analysis from either weak- or strong- 
coupling limits, it is natural to try to reach these points 
as a result of the flow away from the c = 3 critical re- 
gion defined in Sec. IIB.3 [Eq. (18)]. We thus consider 
the effect of relevant and marginal perturbations on the 
corresponding CFT originating from the different inter- 
actions in the Hamiltonian (1). 

B. S'{7(4)-symmetric basis 

It is clear that the most convenient way to obtain the 
general continuum limit is by direct fermionization of the 
X-operators which are the generators of the SU{4) alge- 
bra. Before this we perform a canonical transformation 
of the Xj° and Xf° operators for all lattice sites i, 

l.io = -^(A^ + zX 30 )^^ 10 , 

Xf° = ^(X™- l X™)^X?°, (34) 

from which one obtains [c/. Eqs. (6,7)] 

pii = ^xfx^ 

a/3 

A\ = Y, X f X ^ Bj = J2x^°X^ 1 -XfX^ 1 , 

a a 

Et = A\A U Er = B\B U (35) 

where a, (3 = 0, ... ,3 and a = 1, 2, 3. 

The generators of the sw(4) algebra may be represented 
by four Dirac fermions, 

X f ' = \( c Kp4, a ~4^ c Kp), (36) 



subject to the local constraint X)a=o 4 a Ck - a = ^ n 
terms of these variables the transformation 13 is sim- 
ply a particular case of the U{1) canonical transformation 
for Cfc j0 fermions, 

Cfc,o — > ickfi, 4,o — > ~ ic k,o- ( 37 ) 

The fermionized version of the Hamiltonian (5) then 
takes the form 

i,a,b 

+ 92{c\ a Cific\ +1Q c i+ha + ct c i , a ct +la c i+ i i o) 

+ 94(4,a C i,a + 4+l.a C i+ha) (38) 
+ 55(1 - 4a C ^«)( 1 _ C I+l,a C m,a), 

where g x = \{J\\ + J x + K), g 2 = |(J|| - J x + K), 
.93 = \{J\\-J X -K), .94 = (Jl + 2K), 55 = 2K, and 
7=1. Note that we have retained the variable 7 as an 
additional pcrturbative parameter. 

In the low-energy limit the lattice fermions are ex- 
pressed in terms of left- and right-moving fermions 
i>L,ipR at the Fermi points, 

^c„ jQ ~ exp (~ik F x) + ^ R . a exp (ik F x), (39) 

where x — na, ±k F = ±n/ra, a is the lattice constant, 
and r the filling factor, which is related to the group di- 
mension (2 for SU(2), 4 for SU(4j). In the continuum 
limit the different operators may be classified according 
to their scaling dimensions in the vicinity of the confor- 
mally invariant points. 

Perturbation around CFT with c = 3 

For the particular set of values of the interaction pa- 
rameters J\\/K = 1, J x /K = 0, which corresponds to a 
point in the ^-symmetric plane, the Hamiltonian may 
be written as [Eqs. (35), (5)] 

H = KY,[Pf-\{Et + Er) 

i 

+2X°°X™ 1 ] (Jl + 2K) £ X™ (40) 

i 

From the fermionic representation (36) one observes that 
Eq. (40) corresponds to the Hamiltonian (38) with the 
specific coupling values g\ = g 2 = .95/2 = K, 93 = 0, 
7=1 and 34 = Jj_ + 2K . 

In Sec. IIB.3 we have discussed the soluble model which 
has c = 3 critical behavior in the region with finite rung 
and rung-rung interactions. Because this model is equiv- 
alent to the Hamiltonian (40) when the second term van- 
ishes identically (7 = 0), it would appear natural that 
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the critical region be extended by an increase of the pos- 
itive rung-rung interaction. We therefore consider the 
perturbation of this critical model caused by the term 
(E+ + ) . In the continuum limit it is clear that the 7 
term breaks the SU (4) symmetry, 

4,a C i,b4+l,a C i+hb - ^L,a^L,b1pR, a ^R,b + ■ ■ ■ , (41) 

producing a marginally relevant product of currents of 
the S't/(4)-symmetric WZW model for 7 > 0. The RG 
analysis 47 shows directly that this interaction opens an 
exponential gap and drives the system into a dimcr- 
ized state. A similar effect occurs in the spin-orbital 
model with symmetry-breaking perturbations. 36,37 We 
note also that in the study of Ref. 48 this term appears as 
the continuum limit of the generators of the Tcmperley- 
Lieb algebra. 26 In the present analysis the combination 
±(Ef+E7)-X?°X?° 1 is the generator of the Temperley- 
Lieb algebra in the projected 3-state-per-site subspace 
(further details are provided in Ref. 22). 

We have shown that the S'C/(4)-symmetric model (38) 
has an exponential gap, and this is expected to persist 
up to 7 w 1, at which point the model becomes that of 
Eq. (1) with parameters J\\/K = 1, J x /K = 0. Thus in 
the ring-exchange model (1) not only this point is gapped 
but also the region around it in the ./^-direction. We 
note that the rung interaction [54 in Eq. (38)] induces 
both relevant (scaling dimension 3/2) and marginal per- 
turbations away from criticality. 

One may develop a perturbative analysis for the en- 
tire problem by starting from this S'C/(4)-symmetric so- 
lution. The SU(4)k=i CFT allows a decomposition into 
the direct sum of two SU(2)k = 2 CFTs, a fact used in 
the analysis of the SU (4)-symmetric spin-orbital model 
in Ref. 37. The addition of different perturbations to the 
5L r (4)-symmetric solutions then leads to a lowering of the 
symmetry, and may induce a flow to one of the SU(2)k=2 
components. These correspond to the two second-order 
phase transitions with c = 3/2 at weak and strong cou- 
pling. The spin density is expressed in terms of the SU (4) 
primary fields and therefore, in addition to the uniform 
(k = 0) and staggered (k = ir) parts, has an oscillat- 
ing part with momentum k = n/2. This result appears 
plausible for the problem at hand in that a semiclassi- 
cal analysis 11 reveals the presence of a ground state with 
orthogonal spin alignment on neighboring sites. We will 
not dwell further on this approach here. 

C. c = 1 CFT 

We conclude this section by considering perturbations 
around the c = 1 CFT solution given by Eq. (15). Al- 
though the point J± = and line AB on the phase dia- 
gram of Fig. 2 possess additional degeneracies, the con- 
formal dimensions of the model are unaltered. The WZW 
model with c = 1 contains one primary matrix field of di- 
mension (j, j), and the perturbation caused by this field 
is related to dimerization or alternation, i.e. it breaks 




FIG. 10: Cross-section of the phase diagram (Fig. 2) in the 
plane J± = 0. The solid line marks the projection of the 
plane of Z4 symmetry, and the solid circle the line AB of 
massless, incommensurate solutions. The dashed line repre- 
sents the second Babujian-Takhtajan transition in the strong- 
coupling regime, which separates the scalar- and vector- 
chirality phases. The half-plane in which a massless region 
exists around the point J x = — K is represented by the dot- 
ted line. 



explicitly the symmetry of translation by one lattice site. 
This type of interaction is absent in the initial Hamilto- 
nian (1), and therefore the relevant perturbations due to 
this matrix field are disallowed by symmetry. The only 
possible perturbation is then a marginal current-current 
interaction which is present in the 33 term of Eq. (38). 
Perturbation around the exact solution given by Eq. (15) 
requires that all other couplings be set to zero. This cor- 
responds to the ray J\ \ + J x + K = and J± = in the 
parameter space of the model, i.e. to the fermionized 
Hamiltonian 

H = h^2[( C i,a C i,o4+lfl C i+ha + 4,0^4+1,^+1,0) 
i.a 

+ (1 - 2ct oCi>0 ) (1 - 2ct Hi0 c i+ i l0 )] (42) 

+ 5(4,a C »,o4+l,a C m,0 + 4,0 C i,a4+l,0 C i+ha), 

where g = J|| — J x — K , and with the condition J|| + 
J x + K = 0, which follows from Eq. (38). The general 
analysis of Ref. 47 is now applicable: when g = g$ < 
the symmetry-lowering term is marginally relevant and 
generates a gap, while if g$ > the term is marginally 
irrelevant and one expects a massless region to extend 
along this ray, as represented in Fig. 10. The point 53 = 
represents a type of Berezinsky-Kosterlitz-Thouless tran- 
sition, as shown already from the exact solution in Sec. II. 
From the ray equation Jn + J x + K = 0, the coupling g is 
proportional to Ji 1 , and thus for Ji | > one expects mass- 
less behavior, while for J|| < the dynamical generation 
of a Haldane gap occurs. One may further expect that 
the massless behavior persists for a small region of values 
of Jl around this ray in the plane J\ \ + J x + K = 0. For 
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values of the parameters such that the phases develop a 
gap, incommensurate behavior is expected, in the form 
of a maximum in the structure function S(q) at values 
q = 9max where < q max < n, over a broad region in the 
vicinity of the line AB in Fig. 2. 

D. Universality classes 

The preceding subsections have revealed the rich crit- 
ical behaviour of the system under investigation. There 
exist two critical surfaces described by CFTs with c = 
3/2, the weak- and strong-coupling regimes, and one crit- 
ical region in the universality class of a c = 1 CFT. The 
existence of further second-order phase transitions may 
be excluded on the basis of general arguments. Both the 
weak-coupling and the corresponding strong-coupling re- 
gions of the phase diagram may be obtained by consider- 
ing perturbations around the fixed points of two decou- 
pled spin chains. This corresponds to the flow from c = 2 
CFTs to c = 3/2 CFTs: the c = 3/2 theories represents 
unstable fixed points, and a small change of the couplings 
in the Hamiltonian leads to further flow to the stable 
c = 1 CFT. For the general model (1) with SU (2) symme- 
try, the latter CFT has the lowest possible central charge. 
These considerations are consistent with the Zamolod- 
chikov theorem. As discussed in Sec. IVB, perturbations 
around the S'[/(4)-symmetric point reduce the symmetry 
and open a gap, thus excluding the possibility of c = 3 
criticality. A candidate for this symmetry breaking is 
SU(4) — > SU(3); however, the exactly soluble critical 
model with c = 2, which has SU(3) symmetry, 22 is "too 
far" from the parameter space of the current model, in 
the sense that it has a positive diagonal-diagonal coupling 
Vd d , and addition of interaction terms to this model in- 
duces a flow to a massive phase. We therefore conclude 
that the list of second-order phase transitions presented 
in this section is complete. 



V. SUMMARY AND DISCUSSION 

We now collect all of the information presented in 
the preceding sections. The phase diagram is shown in 
Fig. 2 in 3D form with coordinate axes J\\/K, J x /K, 
and J±/K. We have chosen this convention to high- 
light the phases and phase transitions obtained at strong 
K, as relatively little is known about this regime. The 
parallel dot-dashed lines demarcate the ^-symmetric 
plane, which intersects the horizontal plane on the line 
J x /K = — 1. This line lies parallel to the coordinate axis 
J±/K, and the part —4 < J±/K < 4, represented by 
the thick dashed line between points A and B, forms a 
critical incommensurate region emerging from the exact 
solution obtained in Sec. IIB and discussed in Sec. III. 

The matrix-product ansatz (Sec. Ill) reveals two lines 
of phase transitions, marked by the straight, solid lines in 
the ^-symmetric plane at J±/K = ±4. The line AA' 



at J±_/K = 4 represents a continuum of second-order 
phase transitions from a rung-singlet phase (Fig. 4) to a 
staggered dimer phase (Fig. 6), the existence of which 
has been confirmed by a number of studies. 13,15 ' 16 In 
Fig. 2 we have shown only this line, which is known ex- 
actly, although in fact it constitutes a part of a surface. 
The line BB' at J±/K = —4 denotes a continuum of 
first-order phase transitions from a ferromagnetic phase 
(Fig. 5) to a form of AKLT state (Fig. 9). Both lines 
intersect the incommensurate line at their ends, which 
thus represent two multicritical points (A and B). We 
have deduced (Sec. IV) the presence of another line (sur- 
face) of second-order phase transitions, the presence of 
which is consistent with the Z4 symmetry of the plane 
in which it is in principle known exactly, and the con- 
jectured form of the surface is represented both by the 
shaded triangle in the region close to the origin of Fig. 2 
and by the dashed line in Fig. 10. This phase transi- 
tion separates the strong-coupling vector-chirality phase 
(Fig. 8), dual to the rung-singlet phase, from the scalar- 
chirality phase (Fig. 7), which is dual to the staggered 
dimer phase. While the letters indicating the locations 
of the phases in Fig. 2 refer in a strict sense only to the 
plane of exact Z4 symmetry, the absence of phase tran- 
sitions other than those present on the diagram means 
that the gapped phases obtained for model parameters 
outside the plane are connected continuously to the ex- 
actly known states, and their physical properties evolve 
continuously as the values of the interactions are further 
changed away from the Z4 plane. 

Renormalization-group techniques and perturbative 
approaches applied around the exactly known solutions 
may be used to determine the fixed-point structure of 
the entire phase diagram. The ^-symmetric plane is 
found also to constitute the transition region between 
weak- and strong-coupling regimes, which is centered on 
the star in the phase diagram of Fig. 2. Analysis of the 
relevance of the additional terms in the model perturb- 
ing the c = 3 CFT (Sec. IV) suggest that all points on 
this plane, other than the transition lines of the preced- 
ing paragraph and the incommensurate line AB, have 
gapped excitations. These spin gaps are generated only 
by marginal symmetry-breaking perturbations and are 
therefore small, a result which was confused with gapless 
behavior in the numerical studies of Ref. 14. Within the 
plane, the RG fixed points corresponding to the c = 3/2 
transitions are unstable (Sec. V), and therefore we ex- 
pect a flow to a stable c = 1 theory, for which the ob- 
vious candidate is the incommensurate line bounded by 
the multicritical points A and B. 

There is a certain mathematical and physical similarity 
between the behavior of the 5 = 1/2 ladder with a ring- 
exchange interaction and the S = 1 bilinear-biquadratic 
model described by the Hamiltonian 

#BB=]T(S 4 -S 4+1 ) + a (S 4 .S 4+1 ) 2 . (43) 

i 

For different values of the parameter a this model con- 
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tains a variety of phases and different forms of critical 
behavior. For a = — 1, the Babujian-Takhtajan model, 
it is in the same universality class as the two c = 3/2 
second-order phase-transition lines of the ring-exchange 
model. The Hamiltonian (43) also shows incommensu- 
rate behavior for the parameter region 1/3 < a < 1. 
However, in the bilinear-biquadratic model this is re- 
stricted to a gapped phase, in contrast to the current 
model where the critical incommensurate regime is re- 
lated to the degeneracy of the c = 1 CFT. We com- 
ment that an explicit mapping from the 5 = 1 model of 
Eq. (43) to a S = 1/2 ladder system may be obtained 
by using a composite-spin representation, 18 but that the 
resulting model does not lie in the phase space (Fig. 2) 
of a ladder with ring-exchange interactions (see Sec. I). 

The most accurate and extensive numerical studies 
performed to date for a ladder system with ring exchange 
are contained in Ref. 16. The first quadrant of the cir- 
cular phase diagram presented in this work corresponds 
to the line (J\\ = J±_ > 0, J x = 0) in Fig. 2. We note 
the complete agreement of our analysis with the numer- 
ical investigation for this line: the four phases of Ref. 16 
are those represented schematically in Figs. 4, 6, 7, and 
8, and related under the Z4 transformation as shown in 
Fig. 3. The complete phase diagram allows us not only to 
verify the Babujian-Takhtajan nature of the rung-singlet 
to staggered dimer phase transition but to confirm that 
the scalar-chirality to vector-chirality phase transition 
lies in the same universality class, and in addition to 
specify the location (Jm = K) and nature (Sec. IV) of 
the crossover between staggered dimer and scalar chiral- 
ity phases, which remained unclear from the numerical 
analysis. The second quadrant of Fig. 1 of Ref. 16 is rep- 
resented by the line (J\\ = J± < 0, J x = 0) in Fig. 2. 
Although we have not obtained specific information con- 
cerning this region of the phase diagram, the properties 
of the regions to which it is connected continuously al- 
low us to deduce that the vector-chirality phase should 
be separated from the ferromagnetic phase by first-order 
transitions to a form of AKLT state, which indeed ex- 
hibits the collinear-spin property found by DMRG. Be- 
cause our considerations do not include negative values 
of K, we refrain from comment on the third and fourth 
quadrants of the circular phase diagram. 

Finally, we summarize briefly the relevance of our re- 
sults for the materials and higher-dimensional systems 
mentioned in Sec. I. Experimentally determined values of 
K for cuprate systems, including the ladder compound 
La6CasCu2404i, suggest that the ring-exchange interac- 
tion may in fact be close to the value required to drive 
the rung-singlet phase to a staggered dimer state, and 
thus that staggered dimer correlations may be detectable. 
In two dimensions one expects that a larger value of K 
would be required to find analogs of the more exotic ( i.e. 
chiral) phases of the model, and thus that these would 
most likely be detectable, if at all, as excitations. Films 
of 3 He have been found to offer larger values of K, but 
these appear both in combination with other multiple- 



spin exchange processes and on a more complex lattice 
geometry which may lead to further topological possibili- 
ties for spin configurations. Our results suggest that this 
restricted-geometry system may provide a very rich spec- 
trum of possible phases, but a considerably more specific 
analysis would be required. 

In conclusion, we have analyzed a general model for a S 
= 1/2 ladder with ring-exchange interactions. By inves- 
tigating the exactly soluble points within the parameter 
space we obtain a complete classification of the phases 
and phase transitions in this type of system. Although we 
have considered the minimal model possessing both cyclic 
four-spin interactions and non-trivial exact solutions, we 
find a rich variety of gapped and gapless phases, of first- 
and second order-transitions, and of commensurate and 
incommensurate excitations, all connected by a complex 
renormalization-group flow pattern. The full phase dia- 
gram provides significant additional insight into the types 
of phases and transitions arising in low-dimensional spin 
systems as a consequence of the cooperation and com- 
petition between nearest-neighbor antiferromagnetic ex- 
change interactions and multiple-spin interactions of the 
ring-exchange type. 
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APPENDIX A 

We summarize here the important formulas used in 
Sec. IV, and describe the fermionization procedure for 
the general bilinear- biquadratic spin- 1/2 ladder in the 
limits of two decoupled chains. 

The continuum limit of the SU(2) Heisenberg model 
is described by an SU(2) WZW model at level k = 1. 
This model has one matrix primary operator g a (a — 
0, 1, 2, 3, 4) of scaling dimension (1/4, 1/4). The right and 
left Kac- Moody currents are fields of dimension (1,0) and 
(0, 1). The relationship between these operators and the 
spin-operator density is given in the continuum limit by 
Eq. (25). 

A pair of level- 1 SU(2) WZW models may be repre- 
sented in terms of 4 Ising fields. The operator content of 
these Ising models provides a set of elementary variables 
for constructing the continuum limit of the general ladder 
model. The critical Ising model is described by a c = 1/2 
CFT in the continuum limit. It contains holomorphic 
and antiholomorphic fields, respectively ip(z) and V>(z), 
with conformal dimensions (^,0) and (0,^). The en- 
ergy operator e(z,z) = iil)(z)^)(z) has dimension (5,5), 
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while the order field a(z, z) and disorder field fi(z, z), re- 
lated to the order field by Kramers- Wannier duality, have 
the same dimension (yj, yg). The holomorphic and an- 
tiholomorphic components of energy-momentum tensor 



are T{z) = -^ipdip and T(z) 



\ipdip. 



The representation of the fundamental lcvel-1 SU{2) x 
SU(2) WZW fields is given in terms of 4 Ising models by 
(see e.g. Ref. 44) 



3i 

92 
93 

9'o 
9[ 

9*2 
33 



0-1(720-30-0 + Ml/"2A<3AlO, 
M10203MO — O'lA i 2A i 3O'0, 
O-1M2O-3M0 + Mlf2M3CT0i 
O"lf2M3M0 — /U1M2C3O0, 

O-1CT2CT3O-0 — /UlM2/"3/"0, 
— /il<T2CT3^o - flM2M3CT0, 
-O-1M2O-3MO + Ml^2M3CT0, 
-01CT2/U3/UQ - M1M2C3O0, 



(Al) 



and the expressions for the SU{2) x SU{2) Kac-Moody 
currents are 



Jl = 


5«(^l^0 - 


i>21p3), 


J 2 = 


5«(V>2^0 - 


V ; 3V'l) ) 


^3 = 


^(^3^0 - 


^1^2), 


4 = 


-±i(V>i^o 


+ "02^3), 


^2 = 


-|i(V>2^o 


+ ^3^1). 


^3 = 




+ ^1^2), 



(A2) 



The fields <?j (i=l,2,3) represent the staggered part of 
the spin-density operator, while the Kac-Moody currents 
correspond to its uniform part. 

The operator-product expansions between Ising-model 
fields are 45 



<r(z, z)ct(w, w) 
fj,(z,z)fj,(w,w) 
a(z, z)fi(w, w) 
fjb{z, z)a(w, w) 



\z — w\ 4 

1 



I I 3 _ 

— + — w| 4 e(w;,w;), 



--i|z-u;|«e(iu,tu), (A3) 



7(2 — 10)5^(10) + 7*(z — tD)5^(ty) 

•\/2|z — w| 3 
7*(z — to) 5 ^(ui) + 7(2 — w)i4>(w) 

\/2\z — w\i 



z — w 
1 



+ 2(z- io)T(ro), 
+ 2(z- tS)T(tS), 



and 

ip(z)ip(w) 
, ip(z)'ip(w) 
V>(z)ct(u>, w) 

\/ 'A\Z — W) 2 

V>(zV(w,w) ;_, , (A4) 

l 4>{z)a{w, Hi) 

■ip(z)fi(w, w) 

e(z, z)e(w, w) 

where 7 = exp(i7r/4). 

These relations allow one to compute the OPEs be- 
tween marginal operators defined in Eq. (27), 



z — w 






w) 




1 

W) 2 


7*c(w 


w) 


V2(z- 


w)? 




w) 


V2(z- 


w)i 


ja(w, 


w) 


V2(z- 




1 




z — w Y' 


5 



2 (z)0 2 (w 
Oi(z)0 2 (w 
W>i$i)(z)0i(io 

(lpa4>a)(z)C>2(w 

(V , o^o)(z)0 2 (w 



|Z - 

3 

|z-u;| 4 
-20 2 



|z — w| 
2 

|z — tu| : 



■Oi + 



|z — w\ 2 

-(^2^2) (to) - (^3^3) (to) 



(A5) 



1*- 


w\ 2 


-(ipii>i)(w) 




_ |z " 


w 2 


-(lp2lp2)(w) 


- {ipiipi)(w) 



\z — w\ 
-(V>o^o)(w) 

I 19 ' ' 

z — wr 



a = 1,2,3), 



\z — w\ 2 

Because the scaling dimensions of all operators are 
known around CFT points, one may proceed to the con- 
tinuum limit for the spin Hamiltonian around the limit of 
two decoupled chains described by the WZW model with 
c = 2. There exist two variants of this limit, at weak and 
strong K, and it is convenient to study the continuum 
limit of the general Hamiltonian 

H = Jj,(Sj-Sj + i + T t -T t+ i) 

i 

+^n(Si-Tj +1 + Tj-S i+ i) + JflSjTj 
+Vbi(S < -S i+1 )(T i -T <+1 ) (A6) 
+^ / Rfi(SrTi)(S i+ i-T i+ i) 
+Vb 2 ,(S i -T i+ i)(S i+ i-T i ) 
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where in addition to the leg (Jl), rung (Jr), and diagonal 
(Jd) Hcisenberg interactions we include leg- leg (Vll), 
rung-rung (Vrr) and diagonal-diagonal (Vdd) four-spin 
interactions. The most relevant contribution from the bi- 
quadratic terms arises from the product of the staggered 
parts of the corresponding composite quadratic expres- 
sions, and from the operator product expansion between 
quadratic product of currents with the quadratic prod- 
ucts of the staggered parts, 

Si-Si+i ~ {0 1 +0 2 ) + {-l) x ' a \g a , 
TYT i+1 ~ (O 1 +O 2 ) + (-l) x/a \g' , (A7) 
Sj-Ti ~ i(Oi-0 2 )-(ei+e 2 + e 3 -3eo) 

+(-if/° ■ [g a (J a ' + J a ') + (J a + J a )g' a ] ■■, 



Si-T. 



i -M+l 



- 5(Oi-0 2 ) + (e 1 +e 2 + e 3 -3eo) 
+ (-l) x/Q : [g a (J a ' + J a ') - (J a + J a )g' a ] :, 



and 



: (Sj-Sj+i) :: (Tj-Tj+i) :~ 

3 

A 2 [-4(0 1 + 2 ) + 6^e fc ], 

fe=0 

: (Sj-Tj+i) :: (Sj+i-Tj) :~ 



(A8) 



A 2 [Oi - 50 2 + 3e - 5(ei + e 2 + e 3 ) + 3^e fc ], 

k=0 

■ (Sj-Ti) :: (Sj+i-Ti+i) :~ 

3 

A 2 [Oi - 50 2 - 3e + 5(ei + e 2 + e 3 ) - 3^ e k ], 



where the cut-off-dependent constant A emerges from the 
operator product expansion between different contribu- 
tions. 



These expressions give in terms of Majorana fermions 
the continuum limit of the Hamiltonian, which is sep- 
arated into singlet, triplet, and marginal parts as in 
Eq. (26), with 



m t = Jr-2J d - X 2 (6V ll -2V dd + 2V R r), 

m s = -2,Jr + QJ d -\ 2 {<6V ll + W dd ~W rr ), (A9) 

Ai = -AJ L + J R /2 + Jd + (-W ll + V dd + V rr ), 

A 2 = -4J L - Jr/2 - J D - {-AV LL - 5V DD - Wrr), 



and operators Oi and 2 as defined in Eq. (27) . This rep- 
resentation allows one to study simultaneously the limits 
of weak and strong K (Sec. IV). We note that the signs 
of the four-spin interaction terms in the expression for 
the masses are different from those in Ref. 10 but the 
same as those resulting from the analysis in Ref. 21 of 
the string order parameter. 
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